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PROGRAM DESCRIPTION 


A. IDENTIFICATION 


Program Name 

Programmer's Name 
Date of Issue 

B. GENERAL DEGCRI PlI ^N 


- Croston Exercis Model 
(Univac 1110 Version) 

- R. C. Croston, GE/JSC, Houston 

- August 1<', 1973 


A mathematical model and digital .-nputer simulation of the human 
•ardiovascular system and it3 c^r.trols were developed to simulate 
transient responses to bieycle ergometer exercise. The purpose of 
the model is to provide a method to analyze cardiovascular control 
hypot v ?ses which cannot be easily tested in an animal o human or 
in a spaceflight environment . 


C. USA IE AND RESTRICTIONS 

Machine and Compiler Required - Univac 1110 TSS and Fortran V 

Peripheral Equipment Required - Magnetic Tape 

- Time Sharing Terminal 


D. PARTICULAR DESCRIPTION 

Equations Used and Derivations - See TIP No. 7^1 -MED-2010 
Definition of Terms Used - See Table 1. 

Detailed Description - Equations describing pulsatile blood flows, 
pressures, and volumes for 28 model sections of the uncontrolled 

rascular u are solved, Th<- ilati y system 

model, Figure 1, is co bined with models of the controlling systems 
to simulate transient responses to exercise. Other characteristics of 
the combined model include gravity effects, muscle pu ping, venous 
tone, venous valves, respiratory frequency, and intrathoracic pressure 
effects . 
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E. DESCRIPTION OF INPUT 

(l) Contact V. J. Marks, GE/JSC, Houston to obtain Univac 1110 sign on 
procedures a: <1 authorization :odes. 

(.’) Interactive >perat,ion of the odel will ask the user if he wishes 
to change the standard set of initialization data shown in Table 1. 
If so, the user enters the appropriate index code and new value 
for each variable as directed by the conversational output of the 
program. 


F. DESCRIPTION OF OUTPUT 

(l) Conversational output of the program will ask the user if he wishes 
to change the standard terminal output list of variables. If so, 
the user enters the appropriate index code and new value for each 
• M ! I red •• 1 by » he • nversut i i m out i ui ■ r- ■ v . 

( ) The standard set of output variables consists of time and 8 depen- 

re output at 1 he end of each heart beat. The 
selected standard variables are as follows: 


1 Simulation Time (sec) 

? Heart Rate ( beats /min) 

S Cardiac Output (liters/min) 

U Stroke Volume (liters) 

5 Mean Blood Pressure (mm l'g) 

6 Systolic Blood Pressure (mm Hg) 

7 Diastolic Blood Pressure (mm Hg) 
Respiration Rate ( Breaths/min) 

9 Intrathoraeic Pressure (mm Hg) 

G. INTERNAL CHECKS AND EXITS 


Exit - A normal stop exit occurs when the finish time specified is 
exceeded. 



H. INDEPENDENT SUBROUTINES 


oWIN 

F0N8W 

DELAY 

ALGO 

CVS 

CONTRL 

XIO 

RLKDAT 


Swtt 'h Input 
Function Switch 
Delay 

Integration Algorithm 
Pulsatile Circulatory Model 
Cardiovascular Controlling System:; 
Conversational Input Output 
Initialization Data 


I. COMPUTER PROGRAM LISTING 


(Attached) 
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ELEMENT 

NO. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 
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TABLE 1 
DEFINITIONS 


CONSTANT OR 


MNEMONIC 

DEFINITION 

INITIAL VALUE 

UNITS 

X(l) 

Stressed Vol . , RT. Atrium 

71.88 

ml 

X ( 2 ) 

M II 

, RT. Ventricle 

158.2 

11 

X(3) 

•I II 

, Left Atrium 

29.62 

ii 

X(4) 

•I II 

, Left Ventricle 

147.6 

11 

X(5) 

II II 

, Pulmonary Arteries 

12.08 

11 

X( 6 ) 

•1 II 

, " Capillaries 

16.51 

ii 

X(7) 

♦ •• 

, " Venules 

24.97 

11 

X ( 8 ) 

Total Vol. 

, Aortic Arch 

165.7 

ii 

X(9) 

Inertance Integral 

0.0 

• 

X(10) 

Integral of 

Aortic Arch 

0.0 

irmHg-sec 


Pres3ur«/Brat 



X(ll) 

Inertance Integral 

0.0 

m 

X{ 12) 

Total Vol., 

Thoracic Aorta 

202.6 

ml 

X(13) 

Integral of Carotid 

0.0 

mmHg-sec 


Pressure/Beat 



X( lU) 

Total Vol . , 

Abdominal Aorta 

112.5 

ml 

X(15) 

Stressed Vol., Common Iliac Arteries 

100.3 

II 

X(16) 

II II 

, Legs Small Arteries 

101.2 

•1 

X(17) 

II II 

, Legs Arterioles 

109.3 

II 

x(i 8 ) 

II II 

, Legs Venules 

81.79 

II 

X(19) 

II II 

, Legs Small Veins 

61.72 

II 

X(20) 

II II 

, Femoral Veins 

4.932 

•1 

X(21) 

Total Vol., 

Abdominal Vena Cava 

597.7 

II 

X(22) 

Total Vol., 

Thoracic Vena Cava 

348.6 

II 

X(23) 

Total Vol., 

Superior Vena Cava 

2.587 

II 

X(24) 

Stressed Vol., Lower Carotid Arteries 

25.37 

II 

X(25) 

II II 

, Upper Carotid Arteries 

25.35 

II 

X(26) 

II II 

, Head Small Veins 

0.0 

II 

X(27) 

II II 

, Jugular Veins 

0.0 

II 

X(2 8 ) 

II II 

, Superior Mesenteric 

264.4 

II 



Veins 



X(29) 

It II 

, Inferior Mesenteric 

32.62 

II 



Veins 



X(30) 

II II 

, Portal Veins 

127.3 

II 

X( 31) 

II II 

, Renal Arteries 

24.31 

II 

X(32) 

II II 

, Renal Veins 

82.68 

II 

X( 33) 

Integral of 

Left Vent. Flow/Beat 

0.0 

II 

X(3M 

Inertance Integral 

0.0 

- 

X(35) 

II 

II 

0.0 

- 

X( 36) 

II 

II 

0.0 

- 

X( 37) 

Integral of 

Pulmonary Arterial Pressure 

0.0 

mmHg-sec 

X(38) 

II II 

Arterioles " 

0.0 

•1 

X( 39) 

Not Used 


• 


X(4o) 

II 


- 
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ELEMENT CONSTANT OR 


NO. 

MNEMONIC 

DEFINITION 

INITIAL VAL'JE 

UNITS 

Ul 

X(*4l) 

(ACCMET) Metabolites 

0.0 


U2 

X(*42) 

Not Used 

• 

• 

U3 

X(**3) 

(XN4) Neural Signal 

0.0 

- 

UU 

X( *4*4 ) 

(DA) Alactic Debt 

0.0 

Liters 0 2 

*45 

X(*45) 

(DL) Lactic Debt 

0.0 

Liters Op 

46 

X( 46) 

(AVDA) AVD Lag 

0.0 

C 

*47 

X(*47) 

Not Used 

. 

- 

*48 

X(U8) 

(D0 2 ) 0 Deficit 
(XN3) Neural 

• 0.0 

Liters 0 o 

**9 

X(*49) 

0.0 

2 

50 

X(50) 

Not Used 

- 

- 

51 - 100 


Time Derivatives of the 

above state variables 



101 

QM 

Flow 

from 

RT. Atrium 

Computed 

ml/sec 

102 

QRV 

tt 

99 

RT. Ventricle 

Variable 

99 

103 

QIA 

ft 

99 

Left Atrium 

99 

99 

10*4 

QLV 

19 

99 

" Ventricle 

99 

99 

105 

QJ>A 

99 

99 

Pulmonary Arteries 

99 

99 

106 

<*PC 

99 

99 

Pulmonary Arterioles 

99 

99 

107 

QPV 

•9 

99 

" Venules 

99 

99 

108 

QAA 

99 

99 

Aortic Arch 

99 

99 

109 


Not Used 




110 


99 

99 




111 


Not Used 


Computed 

ml/sec 

112 

QLTA 

Flow 

from 

Thoracic Aorta 

Variable 

99 

113 


Not Used 




11*4 

QLABA 

Flow 

from 

Abdominal Aorta 

99 

99 

115 

QCIIL 

Flow 

from 

Common Iliac Arteries 

99 

99 

116 

QLGSA 

99 

99 

Leg Small Arteries 

99 

99 

117 


Not Used 




118 

(iLOCAP 

Flow 

from 

Leg Arterioles 

99 

99 

119 

QLOVE 

99 

99 

Leg Venules 

99 

99 

l20 

QLGSV 

99 

99 

Leg Small Veins 

99 

99 

121 

QFEV 

99 

99 

Femoral Veins 

99 

99 

V2 

QABVC 

99 

99 

Abdominal Vena Cava 

99 

99 

12' 

QTHVC 

99 

99 

Thoracic 

99 

99 

124 

QSPVC 

99 

9* 

Superior Vena Cava 

99 

99 

125 

QLOC 

Flow 

to Lower Carotid Arteries 

99 

99 

126 

QUPC 

99 

" Upper 

99 

99 

127 

QHCAP 

99 

from Upper Carotid Arteries 

99 

99 

128 

CiHSV 

99 

from 

Head Small Veins 

99 

99 

129 

QfJV 

99 

from Jugular Veins 

99 

99 

130 

QC0R 

Coronary Blood Flow 

99 

99 


iw l\J 



ELEMENT 



CONSTANT OR 

7 

NO. 

MNEMONIC 

DEFINITION 

INITIAL VALUE 

UNITS 

131 

QCSMA 

Flow to Superior Mesenteric Veins 

Computed 

ml/sec 

132 

QIMA 

" Inferior Mesenteric Veins 

Variable 

tt 

133 

3C3MV 

PI ow from Superior Mesenteric Veins 

M 

M 

134 

QP0V 

" Portal Veins 

tt 

tt 

135 

QTMV 

" Inferior Mesenteric Veins 

M 

tt 

136 

tJvENA 

Flow to Renal Arteries 

•t 

tt 

137 

QRALE 

Flow from Renal Arteries 

tt 

tt 

138 

QRENV 

" Renal Vein 

tt 

tt 

139 

QRET 

Flow to RT. Atrium 

tt 

tt 

lUo 

QD(1) 

Not Used 



l4l 


Not Used 



142 


Not Used 



143 


Not Used 



144 


Not Used 



1 U 5 


Not Used 



146 


Not U 3 ed 



147 


Not Used 



148 


Not Used 



149 

QP(10) 

Not Used 



150 

Q3 KB 

Flow Through Skeleton, Bone Marrow, 

Computed 

ml/sec 



and Fat 

Variable 

tt 

151 

CRA 

Compliance, Right Atrium 

tt 

ml/mm Hg 

152 

CRV 

, Right Ventricle 

tt 

tt 

153 

CLA 

, Left Atrium 

tt 

It 

154 

CLV 

, Left Ventricle 

tt 

tt 

155 

CPA 

Compliance, Pulmonary Arteries 

1.2 

ml/rnmHg 

156 

CPC 

Arterioles 

1.7 

tt 

157 

CPV 

" " Venules 

5.3 

tt 

158 


Temporary Storage 



159 


M tt 



160 


H tt 



161 


It It 



162 


It tt 



163 


•t M 



164 


tt tt 



165 

CCILL 

Compliance, Common Iliac Arteries 

0.8 

tt 

166 

CLGSA 

, Leg Small Arteries 

0.8 

tt 

167 

CLGAR 

, Leg Arterioles 

0.3 

It 

168 

CL3VE 

, Leg Venules 

Computed 

tt 




Variable 


169 

CLGSV 

, Leg Small Veins 

tt 

tt 

170 

CFEV 

, Leg Femoral Veins 

0.2 

tt 

171 


Temporary Storage 



172 


M It 



173 


M tt 



174 

CL0C 

Compliance, Lower Carotid Arteries 

0.3996 

tt 

175 

CUPC 

It * f _ It M 

, Upper 

0.3996 

tt 

176 

CHSV 

, Head Small Veins 

5.328 


177 

CJV 

, Jugular Veins 

0.9058 

It 

178 

CCSMV 

, Superior Mesenteric Veins 

9.59 

tt 

179 

CIMV 

, Inferior Mesenteric Veins 

1.509 

tt 

180 

CP0V 

, Portal Veins 

6.047 








E.LENENT 



CONSTANT OR 


. 

•TIEMOMIC 

DEFINITION 

INITIAL VALUE 

units 

181 

CRENA 

Compliance, Renal Arteries 

0.222L 

M 

182 

CRENV 

, Renal Veins 

2.517 


183 

crd) 

Not Used 



;8U 


N 



185 


ft 



l9o 


N 



187 


ft 



188 


M 



189 


H 



190 


•f 



191 


•• 



192 


tf 



193 


ff 



19 1 * 


II 



195 


II 



198 


II 



197 


II 



198 


II 



199 


•1 



200 

CD(l8) 

II 
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"ALCUIATED VARIABLES FOR EACH VOLUME COMPARTMEfT’ 


BUWSNT 
. NO, 

PRESSURE 
(mm H< 




PRESSURE DUE 
TO GRAVITY 


EXTERNAL 

PRESSURE 

VOLUME 

COMPARTMENT 

201 

PRA 

321 

VQ) 

421 

PG(1] 


455 

PEXT(l) 

RT. Atrium 

202 

PRV 

3 22 

V(2) 

422 

•a 


456 

PEXT(2) 

RT. Ventricle 

203 

PIA 

323 

V( 3) 

423 

PG(3 


457 

PEXT(3) 

Left Atrium 

204 

PLV 

324 

V(4) 

424 

P0(4] 


458 

PEXT(4) 

Left Ventricle 

205 

PPA 

325 

V(5) 

425 

P0(5) 


459 

HDCT( ) 

Pulmonary Ar- 

206 

PPC 

326 

V(6) 

426 

PG(6) 

460 

PEXT(6) 

teries 

Pulmonary Ar- 

207 

PPV 

327 

V( 7) 

427 

P0(7) 

461 

PEXT(7) 

terioles 
Pulmonary Veins 

206 

PAA 

328 

V(8) 

428 

P0(8) 

462 

PEXT( 8) 

Aortic Arch 

209 

Temporary Storage 

429 

PG(9) 

463 

PEXT(9) 


210 

211 

212 

•t 

M 

PLTA 

•• 

•1 

332 

v( 12 ) 

430 

1*31 

432 

P0(10) 

w(u) 
po( 12 ) 

464 

465 

466 

PEXT( 10) 

rEXTlipi 

Thoracic Aorta 

213 

214 

Temporary Storage 
PLABA 33U V(l4) 

512 


!! 

467 

468 

PEXTf 13) 
pext(i4) 

Abdominal Aorta 

215 

PC ILL 

335 

v( 15 ) 

435 

P0(15) 

469 

pext( 15 ) 

Common Iliac 
Artery 

Leg Small Art- 

216 

PLGSA 

336 

v( 16 ) 

436 

P 0 ( 16 ) 

470 

PEXT(l6) 

217 

PI/1AR 

337 

v( 17) 

437 

PG(17) 

471 

PEXT( 17) 

eries 

Leg Arterioles 

218 

PLGVE 

338 

v( 18 ) 

438 

PG(l8) 

472 

pext( 18 ) 

Leg Veins 

219 

PLGSV 

339 

V(19) 

439 

PG(19) 

473 

PEXT( 19) 

Leg Small Veins 

220 

PFEV 

3U0 

V(20) 

44o 

rc( 2 o) 

474 

PEXT(20) 

Femoral Veins 

221 

PABVC 

3Ul 

V(21 ) 

441 

PG(2l) 

475 

PEXT(2l) 

Abdominal Vena 

222 

PTHVC 

342 

V(22) 

442 

PG(22) 

476 

PEXT(2?) 

Cava 

Thoracic Vena 

223 

PSFVC 

343 

V(23) 

443 

FG(23) 

477 

PEXT(23) 

Cava 

Superior Vena 

224 

PL0C 

344 

V(24) 

444 

PG(24) 

478 

PEXT(24) 

Cava 

Lower Carotid 

225 

PUPC 

345 

V(25) 

445 

P0(25) 

479 

PEXT(25) 

Arteries 
Upper Carotid 

226 

PHSV 

346 

V(26) 

446 

PG(2o) 

480 

PEXT( 2b) 

Arteries 

Head Small Veins 

227 

PJV 

347 

V( 27) 

1*47 

PG(27) 

48l 

PEXT(27) 

Jugular Veins 

228 

FCSMV 

348 

V(28) 

448 

PG(28) 

482 

PEXT(28) 

Superior Mesen- 

229 

PIMV 

349 

V(29) 

449 

PG(29) 

483 

TEXT(29) 

teric Veins 
Inferior Mesen- 

230 

PP0V 

350 

v( 30 ) 

450 

PG(30) 

481* 

PEXT( 30) 

teric Veins 
Portal Veins 
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BLBUafT 

PRESSURE 

TOTAL VOUWE PRE3SURK 

DUE EXTERNAL 

VOLtME 

. 

(m Mg) 

(ml) TO (SRAVITY PRESSURE 

COMPARTMKlfT 

231 

nuniA 

351 V( 31) **51 P0(3l) 

V' 1 PEXT( 31 ) 

Renal Arteries 

232 

PRENV 

352 V( 32) 1*52 PO(32) 

4(V PKXT( 32) 

Renal Veins 

233 

PD(l) Mean Pulmonary Arteries 



21U 

PD(2) Wedge Pressure 1*53 - 454 Not 

Used 


235 

Not Used 353 - 368 Not Used 



236 

M N 

i6‘ > - 7(1**') Blood Volui •- Command i .0 1 


23 

w I! 

370 - V(50) Total Hl<xxi Volume 

- 5325.0 ml 


236 

H H 




239 

M H 




2U0 

•I N 




21*1 

M H 




21*2 

t« *• 




21*3 

M II 




21*1* 

M »• 




2^5 

M II 




2 !*»', 

•I II 




21*7 

H M 




246 

•I H 




2U9 

PM M' an 

Arterial Pressure 

90.0 

mmHg 

250 

1MC Mean 

Carotid Arterial Pr' ssur' 

90.0 

mmHg 

ELEMENT 



CONSTANT OR 



MNEMONIC 

DEFINITION 

INITIAL VALUE 

UNITS 

251 

RRA 

RT. Atrium Valve Resistance 

0.007508 

mmHg/ml/sec 

252 

RRV 

RT. Ventricle Valve Resistance 0.007508 

M 

253 

RMV 

lie ft Atrium Valve 

0.007508 

N 

254 

RAV 

Left Ventricle Valve 

0.0066 ll* 

II 

255 

RPA 

Pulmonary Arterioles 

0.01602 

•I 

256 

RPC 

Pulmonary Capillaries 

0.05255 

H 

257 

RFV 

Pulmonary Venules 

0.01502 

N 

258 


Not Used 



259 


Not Used 



260 


Not Used 



261 

RLi'A 

Thoracic Aorta 

0.021*00 

•I 

262 


Not Used 

0.0 


263 

RLAHA 

Abdominal Aorta 

0 . 02^*00 

II 

261* 

RCII.L 

Common Iliac Arteries 

0.021*00 

•f 

265 

RLT.SA 

Leg Small Arteries 

0.03003 

•I 

266 

RWAR 

Leg Arterioles 

4.505 

»» 

267 

RLOCAP 

Leg Capillaries 

0.1*505 

II 

268 

RLSVE 

Log Venules 

0.07508 

II 

269 

Rl-T.SV 

Leg Small Veins 

0.07508 

II 

270 

RFEV 

Femoral Veins 

0.02102 

II 
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ELEMENT 



CONSTANT OR 


NO. 

MNEMONIC 

DEFINITION 

INITIAL VALUE 

UNITS 

271 

RAB7C 

Abdominal Vona Cava 

0.007380 

mmHg/ml/sec 

272 

RTHVC 

Thoracic Vena Cava 

0.007508 

II 

273 

RSPVC 

Superior Vona Cava 

0.01502 

N 

27U 

RL0C 

Lower Carotid Arteries 

0.02252 

H 

275 

RUPC 

Upper Carotid Arteries 

0.03378 

M 

276 

RHCAP 

Head Capillaries 

3. **31 

•• 

277 

RH3V 

Head Small Veins 

0.375* 4 

N 

278 

RJV 

Jugular Veins 

O.OOU 302 

II 

279 

RO0R 

Coronary 

15.390 

II 

200 

RCSMA 

Superior Mesenteric Arteries 

2.350 

II 

281 

RIMA 

Inferior MeBenterie Arteriea 

3**.534 r > 

H 

282 

RCSMV 

Superior Mesenteric Veins 

0.2252 

N 

283 

RPOV 

Portal Veins 

0.5255 

II 

28 L 

HIMV 

Inferior Mesenteric Veins 

0.3003 

•I 

285 

HRKNA 

Renal Arteries 

0.01502 

II 

206 

HR ALE 

Renal Arterioles 

0.»450l*5 

•1 

287 

RREFF 

Efferent Arterioles 

2.7L1* 

N 

288 

RRENV 

Renal Veins 

O. 6 L 9 U 

1* 

289 

RD(1) 

Not Used 



290 


•• II 



291 


Not Used 



292 


•I *» 



293 


M •• 



29 1 * 


•I H 



295 


M M 



296 


•I II 



297 


•• »» 



298 


M If 



299 

RD(ll) 

•• It 



300 

RSKB 

Skeleton and Fat 

5.150 


301 

FLPA 

. Inertance, Pulmonary Arteries 

. 

mi'iHg/ml/sec c 

302 

F1.AA 

, Aortic Arch 

0.002881 

II 

303 


Not Used 



30L 


•• M 



305 


II II 



306 

FLLTA 

Inertance, Thoracic Aorta 

0.008538 

•1 

307 


Not Used 



308 

FLLARA 

Inertance, Abdominal Aorta 

0.008977 

II 

300 

FLCILL 

, Iliac Arteries 

0.00626 

•1 


310 - 320 


Not Used 
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UNS 1? Kf.SKD VOLUME." 


UM- .. 

HO. 

MNEMONIC 

271 

vu(i) 


372 

VU(2) 


373 

vu(3) 


37 L 

VU(l/) 


375 

VU(5) 


376 

VU(6) 


377 

VU(7) 


378 

VU(R) 


379 

VU(9) 


380 

VU(lO) 

381 

VU(ll) 

382 

VU(12) 

383 

W(13) 

384 

VU(ll4) 

385 

VU(15] 

I 

386 

VU(l6) 

1 

387 

VU(17) 

388 

VU(18) 

389 

VU(l9) 

390 

VU(20) 

391 

VU(2a; 

392 

VM(22) 

393 

VU(23) 

39^ 

VU(2*4 ) 

395 

VU(25) 

396 

VU(26j 

1 

397 

VU(27j 

1 

398 

VU128 

1 

399 

VU(29) 

1 

*400 

W( 3 o] 

1 

*401 

'n)( 31 ! 

> 

*402 

VU(32J 

1 

*403 - 

*420 


*487 

e(i ) 


*488 

E(2) 


*.89 

E(3) 


L90 

E(*4) 


*.91 - 

500 



DEFINITION' 


CONSTANT OP 

INITIAL VALUE ifNITS 


Rt. Atrium 

30. L 

Rt. Ventricle 

0.0 

Left Atrium 

30.0 

Left Ventricle 

0.0 

Pulmonary Arteries 

. 

Pulmonary Arterioles 

15.0 

Pulmonary Veins 

*400.0 

Not Used 



•t •• 


ml 


M It 

N M 

H •• 

Common Iliac Arteries 
Leg Small Arteries 
Leg Arterioles 
Leg Veins 
Leg Small Veins 
Femoral Veins 

Not Used 
•• ♦« 

II •• 

Low- r Carotid Arteries 
Upper Carotid Arteries 
Head Small Veins 
Jugular Veins 
Superior Mesenteric Veins 
Inferior Mesenteric V' ins 
Portal Veins 

Renal Arteries 
Renal Veins 
Hot Used 

Right Atrial Elastanee 
Right Ventricle Elastanee 
Left Atrial Elastanee 
Left Ventricl Elastanee 
Temporary Storage 


• 

3 . 

30.0 

100.0 

68.0 

*40. 0 


50.0 

50.0 

500.0 

28.0 

*451.0 

. 

375.0 

. 

. 


Computed Variable mmHg/ml 

II II I* 


If II If 
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UNSTRESSED VOLlflgS 


LEML'NT 

U . 

MNEMONIC 

DEFINITION 


CONSTANT OP 
INITIAL VALUE 

501 

z(i) 


Length of Vascular Segment 

0.0 

502 

Z(2) 


It 


0.0 

503 

7 1 

3) 


H 


0.0 

5C»* 

z( 

U 


H 


0.0 

505 

Z( 

5) 


•• 


0.0 

506 

Z| 

6) 


•1 


0.0 

507 

Z 

7) 


It 


o.t 

506 

z( 

8) 


II 


-3.0 

539 

Z( 

9) 


II 


0.0 

510 

z( 

10) 

II 


0.0 

5x1 

Z( 

11) 

II 


0.0 

512 

Z< 

12 


II 


20.0 

513 

z 

13) 


II 


0.0 

51** 

Zl 

1*0 


II 


16.0 

515 

z( 

15, 


II 


0.0 

516 

Z| 

16) 


•I 


0.0 

517 

z( 

17) 


II 


0.0 

518 

Zl 

18) 


II 


0.0 

519 

z< 

19) 


II 


0.0 

520 

zi 

20) 

II 


6.0 

521 

z(2l) 

II 


16.0 

522 

Z(22) 

It 


10.0 

523 

Z( 

23) 

It 


-3.0 

521* 

z( 

2*0 

II 


-1*4.0 

525 

ZI 

25) 

II 


0.0 

526 

z 

26) 

II 


0.0 

527 

ZI 

27) 

II 


-]*4.0 

528 

z( 

28) 

II 


0.0 

529 

z(29) 

II 


0.0 

530 

Z(30) 

♦1 


0.0 

531 

Z{ 

31) 

II 


0.0 

532 

Z{ 

32) 

II 


0.0 

533 

ZI 

33) 

II 


0.0 

53** 

ZI 

3*» ) 

II 


0.0 

535 

r» j 
U 1 

35, 


II 


0.0 

536 

Z 

36) 


II 


0.0 

537 

ZI 

37! 


II 


0.0 

538 

ZI 

38) 


II 


0.0 

539 

Z 

39 


It 


0.0 

5**0 

z(*4o: 


II 


0.0 

5*4 1 

Wh 

*(l) 


Tim"(Sec), Work(KFM/Min) Stepr, 

0.0 


w 

2 


It 

II 

0.0 

5*43 

WK(3) 

II 

II 

8.0 

5*4*4 

WK(U) 

II 

II 

600.0 

5*‘5 

n 

5) 

It 

II 

320.0 

5*x6 

w 

:(6) 

It 

II 

0.0 

5*47 

WK(7) 

II 

II 

0.0 

5>*8 

w 

:(8 


II 

II 

0.0 

5** 9 

w 

9 


II 

It 

0.0 

550 

WK 

(10) 

II 

II 

0.0 


UNITS 

cm 


&6C. 
KFM/Min 
Sec . 

KTO/Min. 

Sec. 

KTW/Min. 

Sec . 

KPM/Min. 
Sec . 

KTO /Min. 



UNSTTtE8SED VOLUMES 


ELEMENT 



CONSTANT OR 


a . 

MNEMONIC 

DEFINITION 

INITIAL VALUE 

UNI'13 

55J 

WK(ll) 

Timo(Sec), Work(KlW/Min) Steps 

0.0 

3ec. 

55? 

WK(12) 

•« M 

0.0 

KPM/min 

553 

WK(13) 

•t M 

0.0 

Sec. 

55^ 

WK(lL ) 

•• ft 

O.d 

KPM/min. 

555 

wk(15) 

•t M 

0.0 

Sec. 

556 

WK(l6) 

M M 

0.0 

KTM/min. 

557 

wean 

tt tt 

0.0 

See. 

55« 

WX(l8) 

•• H 

0.0 

KPM/min. 

559 

VK(19) 

•• •• 

0.0 

Sec. 

560 

WK(20) 

Finish Time 

320.0 

Sec. 

56l 

HR 

Heart Rate 

72.0 

Beat r. /Min 

562 

SV 

Stroke Volume 

0.09 

Liters 

563 

CO 

Cardiac Output 

6.7 

Liters/Min 

56L 

RT 

Total Peripheral Resistance 

Calculated 

mmHg/L/Min 

565 

PEX 

Exercise Boolean (Floating) 

0.0 


566 

W 

Work Rate 

0.0 

KJW/Min 

567 

PSYS 

Systolic Blood Presstve 

0.0 

mmHg 

568 

PDYS 

Diastolic Blood Presture 

0.0 

ft 

569 

FREQ 

Respiratory Frequency 

8.3 

Beats/Mln 

570 

V02D0T 

Oxygen Uptake 

Calculated 

Liters Op/Mi n 

571 

AVD 

Arterio-Venous O 2 Difference 

0.055 

Liters 0g 





Liters Elocd 

572 

PUR 

I ntra- Abdominal Pressure 

Calculated 

mmHg 

573 

PITH 

Intre-Thoracic Pressure 

II 

It 

57I4 

PMP 

L<g Muscl' Pump During Exercise 

ft 

It 

575 

THETA 

Body Angle Relative to riu/ izontal 90. 0 

Degrees 

576 

SF 

Constriction Strength Factor 

1.0 

- 

577 

TTOT 

Heart Period 

0.833 

Sec. 

578 

TAR 

Period of Systole 

0.19 

Sec. 

579 

TV? 

Period of Diastole 

C .36 

Sec. 

580 

Cl 

Gain Constant 

I 16.0 

- 

581 

C2 

Gain Constant. 

10.0 

• 

532 

GNEW 

Gain Constant 

- 0.005 

- 

5> * 

PEX IN 

Pressure Set Point 

88.0 

mmHg 

58 '. 

TR 

Respiratory Period 

Calculated 

sec . 

585-600 


Not Used 





C JO a *F*CtS **npri A I T M OATA IN|T|*Ll7*T|oN */ 7/71 
COMMON/STATr/IISOI , KDOT I SO » 

2/«5TATf/ORA*QPV.OL*.9lV»$PA»QPC»QPV»OAA.QARC,Q^»A.QUTA»QLTA»©tlABAt 

301 ABA ,0C ILL ,9l6SA .OLGA* ,OLGC AP .OlC-VP . «L 6 S V . Of C V . 0 * B V C 18 THVC ,OS*VC, 
14 0 1 0 C » Q H p f , fj m C » p iC'M^V . 0 J V , 8 C 0 • , 0 f S m » , Q I*** iQCSMtf , Q P 0 V • Q I M V • 
SU»tM,QBUriU»FNV,'3*f T *QP< IP) . CSKP 

6 / 5 T»TE/CRA .C*V .Cl* .CIV.CPA.CPC .CPV.C *A.C»HC .Cl AA .CUTA.ClTA.CU***. 

7Cl*P*.CflLl,Cl6S*.ClG* B .ClGVr,ClCSV,CFCV,C*HvC.CTMVC.C5^VC. 

*C10C.CUPC.CM5V.CJV,CCS**V,CIHV,CP0W. 

*C «tN* ,C»CNV ,Cf) I I P ) 

*/5T 1 T[/M»,*»v,ri»,PlV.PB».P*Ci*7v.PM.P<llC.*l»».' , UT».PU*.rU»*», 

B F i. A n A .PC ILL . p l G S A » P l G A P * PL G V C i^LtSV » p F l V * P A 8 V C » P Tm VC * p SP VC * 
CPL0C.PUPC.PM^V.PJV,PCS M V,P|MW,PP0V. 

OPPCNA.PPCNV.POI I*) ,P* » Pp C 
COMMON/^TATT * 

C»PA,OPV,PMV,PAV,RPA,PPC,PPV.»APC|PLAA.PUTA,«tTA.»UA|A. 

r»LAAA,RCILL."L6S».^i.*- 4,, . , L6CAP.PLCVF.RLC,Sv,PrtV.*»*VC» 

CPTHvf»RSPVfi*L0C.*UPC.*HCAP»»HSV#PJV.RC0M l *clMA»PIHA|*C5HV» 

MBPOv.R|“V,HPFN».B»ALC.PPCrF .PRtNV.POl I l ) • R SK B 

l/STATC/rLPM.rLAA.riAPCtFlL^.rLUTA.rLLTA.FLllARAt 

jFLLAPA,rLCTLL,fLCSMA.FLI^A.rLRE*FA,rLOM(B) 

K/sTAtc/vi^ni.vuiso).p<i , 3 M »» p E‘* T ‘3? , *E< ,4 >*SPAci:iin) 

L. ?M0).*ifl?0),MB,SV,C0.BT,PE>l,A,P«;7S. p 0VS,rREQ 

M, V0700T,AVn,P|AB,P|TM. p * , P.TMETA,5r 
M,TTOT,TAS,TVS,CI ,C 2 ,GNf*,PEKlN.TR,DUMPYI I A) 

CALL * I 0 « T ) 

00 |0 I ■ I .32 

| n p< 5< I !«StN(TMfTA/«>7, 2»SB ) *Z I t )•! •0B«98n»/l332. 

CALL CONTRLiT) 

CALL CVS(T) 

c oo 2 n I ■ i , aoo 

C 20 »«» I TF « I np , 30 ) I , * M ) 

C 30 F ORP A T ( S* , I H ( , | 3 . 2 m ) ■ .F I S . 7 ) 

I CALL Al r «OIT) 

IF < T »L T . PK I 20 ) ) GO 70 I 

STOP 

CNO 
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SUBROUTINE c v 5 1 T > 

c c aro i ovascul** rxcRcisr nooel 

C R. C. CBOSTON I 1 f 7 3 SIMPLIFIED VERSION! 

C CONTROLLED SYSTEM 

COMmon/ST ATE / 1 1 SOI • HOOT I so) 

?/STATE/0RA,Q B V,QLA,QLV*QPA,8PC»QPV#QAA,QARC,QLAA,0UTAtQLTAt0UABA, 
3QLABA,QflLLi'9L6SA,OLGA*#QLGCAP,QLr,Vrt«LGSv,QFEV,OABVC»OTMVC,OSPVC, 
*• QL Oc . Q Up C i QHf A P , QMS V » Q J V • Q ( 0 R . Q C S* A ,Q I M A iQCSHV t QPO V » Q I m V » 

B Q B E N A iQbaLF t QR E N V * QRF T * OD 1 10) . QSK R 

6/STaU/CRA,C b V.CLAiCLV»CPAiCPC*CPViCAA,CArc,CLAA.CUTA»CLTAiCUABA, 

7CL»BA,CflU. ,rLG5».CLG* B *CLGVr,CLGSV,crEV,CA8vC.CTMVc.CSPVC» 

BCLOC ,CU B C ,CMSV .C JV ,CCS*V ,C l**v,CPOv i 
f C R E N A iC # ENV » C D I |P ) 

•/STATf/PBA|PBv.PLA,PLViPPAfPPCtPPV»PA»iP»PC»PLAAtPUTA»PLTA»PU*B*i 

BPL*BA t PCILL.PlGSP,PLGP«,PLGVr,PLGSV,PEFV.PABvC.PTMVC,BSPYC. 

C P L 0 C i P U P C ,PHSv *PJV » p C 5 ►* V |P I Mv ,PPOV • 

OPBCnA»PBCNV*PDI IB) »PH|PMC 
COMmon/ST ate / 

EBBA , B R V , B N V |BAV , R P A »BPC » *P V |BABC i BL A A »PUT A * RL T A » Bl! A B A I 
FBL*B*.PClLL.»LGSA.I»LG*P.«»lGCPPiPLGVriRLGSv.PrEV,*»BVC» 

Q R T M v C i R S P V C * * L 0 C t P U P C » R H C A P i R M S V »BJv » R C 0 N i R C S M A t B I M A » B C S M V t 
MBPOViR|MV,BREnA,RRALE i»BEFFt»RENV,BDl I 1 > • PS* B 
I/STATC/FLPA ,FtA A »PlARC iFLLAA^LUTA iFLLTAiFLUA&A* 
jFLLA"A,rLCILL,FLCSMA,FLlMA,FLBCNA,FLOMIA) 

K/STATE/V(Sn)*VU(S0|,PG(3‘O»PC*TI3?),Cl 4 O»SP*CC«|f'l 

l.z< b o).*ki?o>,mr,sv,co.bt,pex,w,psysiPovs,freq 

M, VO?OOT,AVO,P|ABtP|TM,PMB,TMrTA,SF 

N. TT0T,TAS,TVS,CI.C2.6NCB,PEXINiTB,nuMMY»lAI 
0 I Mf NS I ON pb5( | ) , C*p < 3 ? 1 .RSO(SO) • F I NR ( I?) 

EQUIVALENCE ( P R S • p B A I i (CPPI 1 I » C R A ) ,(RSO( 1 ) ,PRA ) , I F | NR ( l / f ”LP A I 
C T IS E L * p SFO T | m £ 

C TT IS A CLOCK FOR ONE BEAT 

TT"T-TSVE 

IF ITT-TTOTI t 00? , I 00 I • I 00 I 
IOr>| TSVC-T 

irST«S p ACF ( I I 
CO" * I 33 I /TTOT* .06 
* ( 3 J I "0 • 0 
PM» i | |0I /TTOT 

x i I n » "0.0 
P M C • X I I 3 I / T T 0 T 
X I I 3 I "0.0 
PD I I I "X I 37 1 /TTOT 
p D ( ? ) • X ( 3 B ) /TTOT 
X ( 37 I "0.0 
* I 3B I "0.0 
SV-TT0T/60. "CO 
RT"PH/CO 

0IFr»-VIS0)"V(R9) . 

X< I B I ■ X I I B I "0 | FF"0.6 
X( I 9 ) ■ X I 1 9 | *0 IFF *0.9 

ps ys^sys 
poTs-ors 
CALL x I 0< T I 
I i n call Contbl « t i 
TEMP-TE“P*n«? 



IF ( TEMP-TI 110,111,111 
III CONTINUE 
STS-0,0 
OTSa | 000. 

TT0Ta«0,/HR » 

TASaO, 10*0.0^*TTOT 
T VSaO, | **0,70*TT0T 

ion? continue 

IFI TT-T AS! 1,2,? 

| SASaSINI 3, ISUiTT/TASl 
f ( I ) -0,0S*0,0S*SAS*SF 
miaO.I ?«0, I *»*5AS*SF 
KSPVC"I?0#t5A5*H0» 1/133?, 

RTHVC" I IO.*SA$a?0*|/|332, 

GO TO 3 
? El I I ■ 0 , OS 
E(3|a0. I? 

RSPVCa.OlSOI* 

RTHVC-.007S07s 
3 TVaTT-O.j 

IF I TV. LT. 0,01 TVaO.O 
|FITV»TVS)N,5,S 
H SVSaS I N I 3, | 1 U*T V/T VS I 
EI2)aO,0|75T.39*SF*SVS 
E I R I - 0.0200M ,Sn*SF»SV5 
GO TO 6 

5 E ( 7 ) *0.01 7S 
E I N | aO ,0? 

6 CONTINUF 

00 II I a I ,<4 
| | CNPI I ) ■ | , /E I I I 

|FI*|NI,LT.0,0)xma0.0 
C COMPUTE VOLUMES 

V I SO I »0 .0 

DO SS 1-1,3? 

VI I l-VUI I ) ♦ X ( I ) 

SS VIS0I-V(50|TV( I ) 

v(%o)-v<50>»vm-vm>-vuj)-v<io) 

C RESPIRATORY PUMPS 

TRSPaTRSP*J-TPS 
I F I TRSP.GT , TR I TR5P"0,0 
T I -TRSP/TR 

P I Th«-? . 67-J9, 70N*TIt5G,M09*T|*TI 
,-S3.R7R*T|*TI*T|tI6.602*T|*T|*T|ttT 
IFlTMETAiGT.AO.IPITHaPITM-1. 

I F I I F S T , E 0 • I IPITH--S. 

P | ARa.P I TH/?, 

DO 71 l-l ,12 

71 PEXTI I JaPlTH 
PEXT ( 22 I »P I TH 
PEXT ( 23 ) -P I TH 
DO 72 I ■ 2 6 , 3 2 

72 PEXT I I I "P I AE 
PEXTI IHIaPJAB 
PEXT I 2 I ) ■ P I A B 

C MUSCLE PUMP 
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TMPaTMP*T-TPS 

TPS-T 

IMTIPiCtil »>TMP«0. 

SHP«51N<?.0,|M4*tmp» 

PMPaRO.aSMP , 

If I SMP . ( T •0#0)PHP»0»0 
IMPf*.lT.|.lf«P*0,0 

in irsi.roi i »php«o.o 

DO HH | m l 4 , | 9 
PC*T II » »PMP 

C # COMPuTC PRESSURES 

00 | ? *1 ■ I » 7 

12 PPS( |)>X< |>/CMPt | ) ♦ PE X T « | » 

00 I) I ■ I S • 20 

13 PRS( | I »x < |)/CmP( I I ♦PE* T I I ) 

DO | H I *2*1 1 32 

IR PPS( tl«Ilt)/CHP( | JpPEKTM ) 

PAAa3.7A4A-0.IOS73*V<8>*O.t)03l7j3#V(fl>»V<e) 

IFIPAA.GT.SYS)SVS»PAA 

IE (PAA.LT, OTS V S ■ P A A 

PLTA*-‘«8.0?B«at42ri8*Vl|7)»0,00n2|S4S*V(l2)«vU2l 
PLAB4«-||.82A*o.On226S*Vll'*)*o.OO»773R*VC|«ti*vcjK| 
PTMVC«-S.sn04*0. I|SH*.I2?)-0,0006S873*VI22I«V(22> 
•«O.OOOOni?34*V(22>*V(22)«V(22) 

PA 6 ^C»- 5 ,‘*R? 6 * 0 . 08 ?‘« 0 e»V< 2 n. 0 . 00 n 33 S , »e»V« 2 J)*vl 21 ) 

.♦0#000000**S07a*V<2|I*VI2I»*V(21) 

PS‘ , VC-*-S.‘«TT’*0.R2S09*V«?3>-0t0‘422R6*V(73»*vi23l 
.♦C, 10063R8S«V(?3)*VC73)«\M23) 

QRAb | PR A»PRV ) /RR A 
C HEART MOOEL 

IEIPRA.LT. PRV)8RA»0.0 
QR Vb X 109) /ELPA 
IE I C°V ,LT .0.0 ) GR V»0 »0 
kootio»»«prv-ppa-rrv*<jrv 

IE(|DOT(OR)»LT.n.O.ANO.OPV.EO.O.O)XDOT(09)a0.n 
OL Aa ( PLA-PLV I / R M V 

If ( pla.lt.plviqla-o.o 

QLVaX I | | )/FL*A 

IE I 0LV.LT.O»O| OLV*0«0 

XOOTIin«PLV-PAA*PGIB)*RAV»QLV 

|EI X DO T ( I I ltLT.0.O.AND.QLV.EQ.0.Q)X00T ( I I | a 0 , n 
C PULMONARY CIRCULATION 

OPAa { PP A-PPC I /RP A 
0 P C ■ ( PPC-PPV I /RPC 
QPVa C PPV-PLA » /RPV 
C ARTERIAL MOOEL 

6 A A a X ( 3RI/ELLTA 

X00TI3R»«PAA-PLTA»PGI12)-RLTa*0AA 
QLTA«XI 3SI/ELLAB \ 

X00T<35»»PLTA-PLAHA*PG(m-RLABA*QLTA 
0LABA-X( 381/fLCILL 

X00T(36)»PLABA-PClLL*P&<IS>-RClLL*QLAb* 

C LEGS 

QClLL-(PCILL-PLGSA)/RLGSA 

QLOSA»(PLGSA-PLGAR)/RLGAR 

QLGCAP-IPLGAR-PLGyEI/PLGCAP 

RL&vr».07S07S 
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IF(8LGVE.LT.C.O)RLGVC»A7 t S67SA7 
QL&VE* I PLGVE-PLGSV ) /RLGvF 

IF I 8 L G V E 1GT*0*0*AN0*FE1«GT iO«0*INn. |FST • C g , | ) r J l G V F ■ I 
RLGSV-.07S07S 

JF(glGSV.LT.O.O)RLGSV»67,SA7SA7 

QLGSV«(PLGSV-PFEV)/RLGSV 

|F(QLGSV.GT»0.0#AND.PEX.GT.0,0tAND.|FST.E8.IIQLGSV»l 

c venous model 

RFC V* »02 I 02 l 

IF(8FEV.LT.0.0>RFEV»67.S67S67 

GFCv«(PfCV-PG< 20 )-PABVC)/RFEV 

8ABVC»CPABVC-PG(2I>-PTMVC)/RA0VC 

0THVC»<PTHVC-PGI22I-PRA)/RTHvC 

Q5PVC«<PSPVC-PG<23)-PRA)/RSPVC 

c head*arms 

8L0C"(PAA*MG«2*n-ML0C)/BL0C 

8UPC«(PL0C--PUPC)/RUPC 

QHCAP»|PUPC-PMSV)/RMCAP 

qh5v»<phsv-pjv)/rhsv 

RJV.,00‘« JOl 

IF(QJV.LT.O.0)RJV-67.S67S67 

QJV«|PJV-PG(27)-PSPVC>/RJV 

C coronary circulation 

QCOR> ( PA A.PRA ) / RC OR 

c CONTINUITY for venous return 

qret»qspvc*othvc*qcor 

c HEPATIC-SPLANCHNIC circulation 

QCSmA»IPLTA-PCSMV»/RCSMA 
OIMAp I PLTA-P I MV I /R IMA 
8CSmv»«PCSMV-PP0V)/RCSMV 
0P0V«(PP0V-PThvC)/RP 0V 
e!MV»(PlMV-PPOV)/R|*V 
C RENAL CIRCULATION 

ORENA»(RL*BA-PRfNA»/RRENA 

QRAlE»(PRENA-PRENV)/(RRALE*RREFF> 

8RENV-IPRENV-PABVCI/RRENV 

C SKELTONiBONE marror,and other 

QSKB«IPLABA.PABVC!/RSKB 
C STATE VARIABLE DERIVATIVES 

XDOT ( I I •8RET-QRA 
XOOT ( 2 ) ■GRA-QRV 
XOOT ( ) J ■OPVOL A 
XOOT ( H ) »0L A-QL V 
XDOT I 5 ) -OR V-QP A 
XOOT I A ) >OPA-QPC 
XDOT ( 7 )»PPC-QPV 
XDOT(6)«8LV-OAA-OCOR-QLOC 
XDOTIIOIaPAA 

XD0T(12)"QAA-8LTA-8CSMA-QIMA 

X DO T ( I R I ■ QL T A»8L*RA-8RENA«8SKB 

XDOT ( IS)«8LABa-8CILL 

XDOT ( I 6 ) «8C ILL-OLGSA 

XDOT ( |7|»0LG5a«8LGCAP 

XDOT (18) >8LGCAP-0LGVE 

XDOT ( 19). 8LGVE-8LGSV 

XDOT ( 20)-8LGSV-8FEV 

X00T(21)«0FEV-8ABVC*QRENV»QSKB 


1 S*8LGvE 


, I S • 8L GS V 



*D0T122)»QABVC^QP0V-QTMVC 
MOOT! 23) aQ JV«QSPVC 
1DOT I 2*» ) -OLOC.QUPC 
I 2S ) •QuPC-QMC *P 
(DOT I 26) «OhC AP-QHTV 
MOOT ) 27 ) •OH5V-Q JV 
* 00 T ( 2D ) •OCS"*-OCS*iV 
*0OT«2^)»0|M*-Q!MV 

*dotoo>«qi*v*qcsmv-qpov 

MOOT ) 3 I ) •ORfN*.OH*LE 

MOOT ( 32 ) •QM*IE«0*ENV 

MOOT t 33 ) ■Ql V 

MOOT It3) aPLOC 

MOOT « 37 ) -PPM 

MOOT ( 3« ) »PPC 

RETURN 

END 
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subroutine COnTRKT) 

COMMON/ ST ATI /I I SO > ,'XOOT t SO > 

2/STATC/QRA,QRv,QIA,QIV,&RA,QRC*9RV»QAA,QARC.01AA,QUTA,QITA,quABA, 

3QlABA,QCIll,QlGSA,QlGAR,QlGCAR.<}l6Vr,UlGSv,QFCV,CABVC,QTMVC,QSRVC, 

••QIOC,9URC,QMCAR,QMSV,QJV,QCOR,QC5mA,Q|MA,qcSmV,QROV,qImv, 

sqrcna , q r a l r , ore n v ,qrft • q r> < tot , w s * r 

A/STATt/CRA,C«V*CLA,CLViCRA,CRC»CRV#CAA,CARc,ciAA,CUTA.CLTA,ci|AbA, 

7CLABA ,CCILL,CU6SA .CLGAR ,CLSVF .CLCSV ,CFCV ,CABVC tCTHVC ,CSRVC , 

• Cl or ,CURC ,CMSv tCJV.CCSMV ,c IMV.CPOV . 

7CRENA tCRFNV.COI IB I 

A/$T A T £ / R R A , R R V » R L A ,PIV » R R A * R R C t R R V • R A A »P»Rf ,RLAA » R U T A » R l T A » P u A 8 A , 
BPlA8A,PClll,ML&SA,PlGAR,R|.GVC,RlGSV,PFtV.RABvr,RTMVC,FSRVC» 

CRIOC ,PUPC ,RMSV ,R JV ,RCSMV ,R IMV ,RPOv , 

ORRCN A , RRCNV.PO I 1*1 ,RH»RHC 
COMMON/ST ATE/ 

FRRA ,RRy ,RMV » R A V ,RRA ,RRC ,RRV * R A H C , R L A A , R U T A , H L T A , R U A B A , 
rRlABA,RClll,BLGSA,NlGAR,RtGCAR,RLGVF|NlGSv,NFEV,AABVCi 
GRTMvC.RSRVC,RlOC.«URC.NHCAR,RMSV,RJV»«COR,RCSMA,R|MA,RCSMV. 
mRROv,R|MV,RRCNA,RRAIE,RREFF , RRE N V , »D ( I I ) i R SK B 

l/STATC/riRA,FlAA,riARC»FllAA,riuTA,FlLTA,FLUABA, 

jFllARA,FieiLl,FLCSMA,FL|MA,FlRCNA,FLDNtB> 

X/STATE/VISOI , V U ( S 0 I , R G I 3 H I »RCXT(32) ,€(**! » S R A C E ( I 0 ) 
l»ZI<«ni.»r«?0 », MR, SV, CO, RT, RE*, R.RSVSiPDYS.FREQ 
M,VO/OOT,AVD,R|AB,RITM,RMP,TMFTA,SF 
N»TT0T,TAS,TVS,C1.C2.GNEA,RE*IN»TR,0uMMYCIa, 

REAL NUM»DI,NUMVO 
EQUIVALENCE 

2 (AcCMET.IIRl I I ,(XNN»X(R3) I , ( D A , x ( R *, I ) , ( 0 1 , X ( R S I I , 

3(002, X ( R B I I , ( X N 3 • 

DIMENSION AVDTSISOI ,V02TSIS0) 

DATA DTS/,001/ 

C SAVE OLD XDOT I R I I 

DO 10 I I ,R* 

10 SPACE I I-3VI«*D0T ( I I 

C *■ WORK RATE FG-M/MlN 

DO 2S |a| , |7,2 
IF ( T-RK (II) 2A , 2S , ?S 
2S »a*K II ♦ | I 
1 4 I F ( R I 27,27,26 

27 REXaO«0 
GO TO 29 

28 REXa | .0 

29 CONTINUE 

C OXYGEN REQUIREMENT FUNCTION V02ROT 

IF(002,lT.O,0)D02an,0. 

V02«DTa.OOOR8SOBlS*#/»2S 
RSRa- 1 , 5*002 
DT l«002 

DT2a ( 2,*002-l ,27S I / I . 15 
DTjar)A*Ol 

DTlN*SR!N(PSR,DT| , 0 T 2 I 
DTaFCNSR(PEX,DT3,OT3,0T|N) 


C 


DAIH-* 1S*(DT-1 ,S)*I .5 


ALACTIC OXYGEN OF B T DA 
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OA|>S*|N(RSRi0O2iOA|H> 

OAO«FCNS*fRC> .0*0.0.0.041 | 

T MFC Nt* I PE I *0.0 .300. .2* I 
ADOT(HRI«IOAO-OAI/T« 

c • LACTIC 0«TGIn OEM PL 

ol ih«.smiot*i .&> 

OL laSMNtPSR.O.O.OLIH) 

OLOaFCNSUMREA.O.O.O.O.DLI » 

T7aFCNS*IPER#0,0.300..2.l 
1 00 T ( RS I ■ ( OLO-OL ) /T 7 

C . aRtcm i al-vcnouS oktgen difference *vn 

NUP9D I •,O3A»0O2 

CALL OCLATIO.O.OS.NUHTOI » AVDTSiNUH»D t » > 

A V 0 • I ( *4 1 *0*OSS 

ROOT ( H4 > • ( NU*«9D- R « HA I I /5 . 

IF(PER.E0.n.0|FLAG»0.0 
IF (FLAG.EQ. I *0160 TO 40 
I F I PC X I 40 • 40 • 4 I 
41 A NF ■ | • 0 

ton«t*?o* 
flag* i . o 

40 IF(T.GT*TONlANF*CtO 

TAN»FCNS«t ANF.3. .34. *3.1 
*DOTIRTI«( I l.00*ANF«RN3>/TAN 
AD0T(H3l»(S.5*P£R-lNH>/4. 

PMHR-2.0 

D m *002*2S«/23, 

IF «P£R,GT.O.OlOTS»OT 
IF (RER.GT.0.0)DHS«OM 
IFIRCR.LT. | . I0H»0«S/0TS*0T 
C ME MON »00 2 / 0 « A 
|F(CMENON.GT.O.S)CHEMON«n.S 
FN*2.*RNH*XN3 
IFIFN.GT.il. )FN> I I . 

IMREI.GT.O. I F NS ■ F N 
|F(PCR.LT.|.lFN«FNl/DTS*OT 
C CONTROLLED RESISTANCES 

RR*LC"(400j*( 1 1020,»SIH0 ,|«Dl/0.T)/|332» 
R5RB*l4i40.*(2|TOO.-A«60,»*OL/0, TO 1/1332. 
RCSHA«l3l30.*IB770 t - , *l30.l*IFN/||,/2.*0H/DMMx/?, 11/1332. 
RMCAR»(*tS70.*(74T0.-A’T5,)*0H/0MMR 1/1332. 
RCOR«(20SOP.-(20SOO.-B770 .)*OM/OMmR)/ 1332, 
RL®CAPa(2AO»»ACCMCT*OAS.OHO.»3RO»*ON/PMx)/|332. 

RL6ARH"! 240. -ACCHET*045.*3RO,.3RO.*OM/OHMx 1/1332. 
RLGarN»«FN*SHO 0»/I I . / l 332 . 

RLGARbS*(00./I 332..RLGARN.RLGARM 

C PRESSURE REFERENCE FUNCTION PR 

PRN«90.*CI »002*C2* ACC-HET 
E3-E7 
E 2 ■ E 0 
EO-EN 

EN*PRN-RM/2.-PHC/2. + RN3»3.*RN<l*FN/2. 

ER*IE3*E2*C0*EN» /Hi 
OOP»O.S33*(ER*GNF.R) 

IF (OOP.LT. O.0)DOP»O.O 
T0T»O. 300*DDP 
MR-40./T0T 
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C CONTROLLED COMPLIANCES 

ERCaPRN-PMC 
IE (C»C.LT.n#0»G0 TO 7 
IE (fRCtGT.NO* )C"C>RO* 

CL&vr»l.TSM(|.0*0 t 0l67ilKCl 
CL^IVO. |R3S« I I .0-0.0167RERC ) 
vu ( i • ) • i on . -E »c 
VU ( |R)»<.e,-E«C*.6ft 
7 CONTINUE 

DELT A-O.R* ( 002-0. RR I 
lE(opl.LT.n.*)DELT»«DO?/6. 

SF b 0,47*ACCMET/R.*DELTA 

C RESPIRATION 

V 0 20 0 T • A VD • C 0 
ERC0«V0200T*P*2R*^.2a 
IE lERE0.6T.in, If RE Q" 30. 

tr»ao./freo 

XOOTIRI !•! V0200T -0.3*1/300. 

IElPEE.ro. O.OIIOOTIHII^-I./JOO. 
IE(ACCMET.LE.0.0.ANO.PEl,EQ*0.0)IO0Tl < (|l |l 0.O 
C 0EV6FN DEFICIT FUNCTION 002 

CALL 0ELAY(O.0.S.VO200T,V02TS,VO2O0,l» 
IOOTIR8)al-VO2OO«PEX*VO2*0T*O. 331/60. 
IF(DQ2.LC.0.0.AND.PFX*EQ,0.O)XDOTf«te>*0,0 
00 31 I»R|,RT 

31 XI I I • X I I 1*0. !•! FOOT! I )♦ SPACE ( |*3t) I 
RETURN 
END 


SUBROUT I NF 0ElAYIFIC.N*X,TS.XOUT.<) 
C NrNO • OF SECS. DELAY 

DIMENSION T S I 1 00 > 

STR0.2 

Mr I F l X I FLOAT I N > /ST I 
IF(K) 10,10,20 
20 X0UT-T5I I I 
DO | I • I ,M 

1 TSC I I • TS I I ♦ I I 
T S I N I ■ X 
RETURN 

10 DO 2 I • » .H 

2 TS I I I «E |C 
RETURN 
END 



subroutine alcoiti 

INTEGRATION ALGOR I T mm 
COMMON /STATE/ I ( SO ) i AOOT < SO I 
DIMENSION IDSl&OI 
IMHI I »| ,7 
7 DO 3 1*1 iNOSV 
3 ADSl | I • I DOT ( I ) 

T • T • M • 

CALL CVS IT) 

00 4 I • I iNOSV 

H « (| >aM/2,R I «00T I t I ♦AOS I | I I ♦« I | ) 
10 RETURN 
I CONTINUE 
NOS V ■ 3 A 
M • • 0 0 2 S 
RETURN 
END 


EUNCT I ON SR I N ( X , A ,M 
IE I A I I ,7,7 

1 SRlNaA 
RETURN 

2 SRlN«B 
RETURN 
END 


E UNC T I ON ECNSR ( * , A ,B ,C ) 

IEI A I I .2,3 

| E C N s R a a 
RETURN 

2 E C N 5 V * B 
return 

3 ECNSR-C 
RETURN 
END 



25 


SUAROUTINF I I 0 I T ) 

COHHON/5TATE/X(a6o> 

OAT* K V , | N | T/ l HV ,0/ 

DATA N2,N3,NA,n6,NA,n7,nP,n9/SA|,S63,SA2i2a9,SA7,SA8,SA9,*7:|/ 

If ( I N | T , 6 T . 0 ) 60 TO 200 
I N I T* I 

ARITE I A i I 0 I 

10 f0AHA»T c *000 rou AISM TO CHANGE 1 N t T I A L I Z C o DATA" I Y / N I * I 
Rf *0 I 5,20) K 
?o for* a n i a 1 1 

If ir.NE.KYl 60 TO AO 
*R I Tf I A, 301 

30 fOAMATI »0f»tt*5C ENTER I NOE * I I - AOO ) , VAlUf, C Rl IIMlIitl'l 
AO R E A 0 I S.iOl I , VAINC A 
50 f 0 Ah a T I I 3 ,r I 2, A I 

If ( I • t T • | .OR, I.6T.AOO) 60 TO AO 
III) ■ VALNf A 
60 TO AO 

60 ARIrC I A, 70) 

70 f ORh A T I • 000 YOU AI5H TO HODlfY THE OUTPU* LIST" l«7M*l 
ACAO I 5,201 K 
If ir.NF.KV) 60 TO 200 

aaitt i a , 6 o i 

• 0 fORMAT MOPLEaSE ENTER A05 I T f 0 A 1 2-7 ) » I Np£ 1 1 I - AOO ) , c*l 112, fAI'l 
90 A E A 0 I 5,|00) IP, | 

100 f OR maT I I 2 , I A I 

If I I A • L T • 2 .OR. I P , 6 T . 9 ) 60 TO 200 
If II.IT.I .OR, I.6T.A00) GO TO 200 
GO TO (90, 102, |03, IOA, 105, 106,107, 108, 109 ) ,:r 
I 02 N 2 ■ | 

60 TO 90 
103 N 3» | 

60 TO 90 
I ON N A • | 

60 TO 90 
105 NS ■ | 

60 TO 90 
I 0 A N6-I 

60 TO 90 
I 07 N 7 • | 

60 TO 90 

108 N8 • | 

GO TO 90 

109 N 9 • | 

60 TO 90 
200 CONTINUE 

If (T.LT.290.) RETURN 

ARITEI A , 2 I 0 I T ,I(N2) ,X!n3) ,I(NA) , X ( N 5 ) , X ( n A ) ,I|N7) , X I N 8 I , X ( n 9 ) 

2 | 0 ' 0 A M A T tf7.2i8F8.3) 

RETURN 

END 
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BLOC* DATA 

cohnon/state/a IIOO) 

<0H*»0N/STATE/B(S01 
CO"*ON/S' ATE'C (SOI 
CO"«ON/ST ATE'D ( SO I 
coM*oN/sT*Tr /f jsoi 
CO«NON/ST*TC/r ( 20 I 
COHhON/STATE/GI 2*01 
C** STATE 

OAT* A/Tt.BB,|SB,?,79,A?,|s7, 6 ,|2.n*,lA.S| ,2«*.97»|AS,7,n # ,P,, 

1 P».207.A.n..||2.S.|OQ.3.|OI,2<IP9.3tB|.79,A| ( 72.M.9p2. 

2 S97,7 i 3mB.*.2.S*7 # 7S.37,2S.3S.0«,0.»2**«,9,37.a2iI27,3, 

3 2 *» . 3 I , A 7 • Afl • A • 0 , , | 0* 0 . • SO* 0 . / 

e*» no* 

DATA B/%Q«p'/ 

c** cohp 

0ATaC'm*O.,I,2,I.7,S.3,7*0.,.B,,p,,3,O.,o.,.?, 

I 3*n,.»)99A l# 399A ( s«32B,,90Sd,9,S9,| f S0S»6,OM7t»?72M,2 # S|7,|«*C»/ 

C»» Pdfs 

Data d/mb*o. .90. »9p./ 

C«* dES 

Data E/3*.O0 7SO*,.OOB*3<«,.0IS02,.nS7SS,.0lsn27,3»O., 

I •02*«iO»,.P2*‘,,02‘‘,«03003 l , *.SOS,,HSoSi,p7SOB,tP7*0»,.02l02i 
2.007 3a ,.022***, #0**S1 3. 1 022S2 , ,033 7* , 3 . «0 | , , 3 7Sd . . 00* 302 , l S . 39 , 2 . 3S 

3 »3‘*,S3**S,,22S2,,S7SS|.3003 i ,0ISP?i, , *S0‘*S,7,7‘«m,.6‘«V*.,0,.0,, 

1 ’*0. |S. IS/ 

C*» | NdT 

OATA f /.OOOTSpp , , 00 ?«B I »3*0.,.0OBs3a,0,,.OPB977,,POA?*,| 1*0./ 
c • • MlSC 

OATa P/**a«0* » S 3 2S . , S32S . , 30 . .0# » 3P . .0. , BS , , i % # , «,pn , #3*0. . 

1 9*0. » S . I 9 N t 3 0 • .30. i I 00. , A B < , N 0 . . 

2 3*0. »SOj .SO. , S 0 0 . .28. • ** S 1 . • *0 » ,37S, • 

3 SO. i ISP. . 16*0. . 3** . 0 i .32*0. ,9*0. . | 0 • 0 . t 

*• 7*0. » - 3 t , 3*0, , 2P. . 0. t I*. . S*n, | a, . I A . . 10* , . 3 , , » 1 9 , , 2*0, , • | d , , | 3*0 , 
S.0.,n.,a,,A0n.,32P,,|M*n,,32p.,72,,.0’ A.7,5*P,,«»3,p., 

* ** # .oSsn l to,ii # , i *33».ia,.3*,96,,in,»-,oos,dH.,i7*o,/ 

END 


I- 10 

1 1 - 20 

2 I • 30 
3»-100 

I 0 I - I so 

isl-l 70 
17 1*200 


2S I -7*0 
2*1-270 
27 l-?*o 
2* I -7»0 
29 1-300 

3(31-320 

32 I - 3 B 0 
3B I - 390 
3 V I - * Ci o 
90 1 -SOO 
So I — S m o 
S H I -S70 
S 7 1 - a 0 (j 



